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The production rate of bottom quarks at hadron colliders exceeds the expectations of 
next-to-leading order perturbative quantum chromodynamics. An additional contribu- 
tion from pair-production of light gluinos, of mass 12 to 16 GeV, with two-body decays 
into bottom quarks and light bottom squarks, yields a differential cross section for bot- 
tom quarks in better agreement with data. The masses of the gluino and bottom squark 
are restricted further by the ratio of like-sign to opposite-sign leptons at hadron colliders. 
Restrictions on this scenario from other data are summarized, and predictions are made 
for other processes such as Upsilon decay into a pair of bottom squarks. 

1. Introduction 

The cross section for bottom-quark production at hadron collider energies exceeds 
the central value of predictions of next-to-leading order OjJLO) perturbative quan- 
tum chromodynamics (QCD) by about a factor of twoBa This longstanding dis- 
crepancy has resisted fully satisfactory resolution within the standard modeltl The 
NLO contributions are large, and it is not excluded that a combination of further 
higher-order effects in production and/or fragmentation may resolve the discrep- 
ancy. However, the disagreement is surprising because the relatively large mass of 
the bottom quark sets a hard scattering scale at which fixed-order perturbative 
QCD computations of other processes are generally successful. The photoproduc- 
tion cross section at DESYls HERAa and and the cross section in photon-photon 
reactions at CERN's LEPu also exceed NLO expectations. The data invite the 
possibility of a contribution from "new physics" n 



2. Supersymmetry Interpretation 

The properties of new particles that can contribute significantly to the bottom 
quark (6) cross section are fairly well circumscribed. To be produced with enough 
cross section the particles must interact strongly and have relatively low mass. They 
must either decay into b quarks or be close imitators of 6's in a variety of channels of 
observation. They must evade constraints based on precise data from measurements 
of Z° decays at CERN's LEP and SLAC's SLC, and from many lower-energy e + e~ 
collider experiments. The minimal supersymmetric standard model (MSSM) is a 
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favorite candidate for physics beyond the standard model. It offers a well-motivated 
theoretical framework and is reasonably well-explored phenomenologically. 

An explanation within the context of the MSSM can satisfy all of the stated 
criteria.cl The existence is assumed of a relatively light color-octet gluino g (mass 
~ 12 to 16 GeV) that decays with 100% branching fraction into a bottom quark b 
and a light color- triplet bottom squark b (mass ~ 2 to 5.5 GeV). The g and the b 
are the spin-1/2 and spin-0 supersymmetric partners of the gluon (g) and bottom 
quark. In this scenario the b is the lightest SUSY particle, and the masses of all 
other SUSY particles are arbitrarily heavy, i.e., of order the electroweak scale or 
greater. (The b either lives long enough to escape from a typical collider detector 
or decays promptly via R-parity violation into a pair of hadronic jets.) Improved 
agreement is obtained with hadron collider rates of bottom-quark production, and 
several predictions are made that can be tested readily with forthcoming data. 

2.1. Differential Cross Section 

The light gluinos are produced in pairs via standard QCD subprocesses, dominantly 
9 + 9 ~~ * 9 + 9 a t Tevatron and Large Hadron Collider (LHC) energies. The g has a 
strong color coupling to 6's and 6's and, as long as its mass satisfies rrig > rrib +mz, 
the g decays promptly to b + b. The magnitude of the b cross section, the shape 
of the b's transverse momentum pxb distribution, and the CDF measurement^] of 
B° — B° mixing are three features of the data that help to establish the preferred 
masses of the g and b. Shown in Fig. 1 is the integrated prb distribution of the 
b quarks that results from g — > b + b, for m§ =14 GeV and m~ h — 3.5 GeV. The 
results are compared with the cross section obtained from next-to-leading order 
(NLO) perturbative QCD and CTEQ4M parton distribution functions (PDF's),! 
with m;, = 4.75 GeV, and a renormalization and factorization scale ii = \J m\ + p\ h . 
SUSY-QCD corrections to bb production are not included as they are not available 
and are generally expected to be smaller than the standard QCD corrections. The 
g-pair cross section is computed from the leading order (LO) matrix element with 
NLO PDF'sf] fj, = ^/t77~ + p\_ , and a two-loop expression for the strong coupling 
a s . To account for NLO effects, this <?-pair cross section is multiplied by 1.9, the 
ratio of inclusive NLO to LO cross sections.i 

A relatively light gluino is necessary in order to obtain a bottom-quark cross 
section comparable in magnitude to the pure QCD component. Values of m~ g ~ 
12 to 16 GeV are chosen because the resulting g decays produce prb spectra that 
are enhanced primarily in the neighborhood of p™™ — rrig where the data show 
the most prominent enhancement above the QCD expectation. Larger values of rrig 
yield too little cross section to be of interest, and smaller values produce more cross 
section than seems tolerated by the ratio of like-sign to opposite-sign leptons from 
b decay, as discussed below. The choice of has an impact on the kinematics of 
the b. After selections on p^ n , large values of reduce the cross section and, in 
addition, lead to shapes of the prb distribution that agree less well with the data. 
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Fig. 1. Bottom-quark cross section in pp collisions at \/S = 1. 
gluino of mass mg = 14 GeV and a bottom squark of mass mg 
is the central value of the NLO QCD prediction. The dotted curve shows the px spectrum of 
the b from the SUSY processes. The solid curve is the sum of the QCD and SUSY components. 
The shaded band represents an uncertainty of roughly ±30% associated with variations of the 
renormalization and factorization scales, the b mass, and the parton densities. The rapidity cut 
on the b's is < !• Data are from Ref. 1. 



After the contributions of the NLO QCD and SUSY components are added, 
the magnitude of the bottom-quark cross section and the shape of the integrated 
distribution are described well. Very ©pod agreement is obtained also with 
data from the UA1 experiment (not shown) .□ The SUSY process produces bottom 
quarks in a four-body final state and thus their momentum correlations are different 
from those of QCD. Angular correlations between muons that arise from decays of 
b's have been measured The angular correlations between 6's in the SUSY case 
are nearly indistinguishable from those of QCD once experimental cuts are applied. 

The energy dependence of the bottom cross section is a potentially important 
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constraint on models in which new physics is invoked to interpret the observed 
excess bottom quark yield. Since the assumed g mass is larger than the mass of 
the 6, the g pair process will turn on more slowly with energy than pure QCD 
production of bb pairs. The new physics contribution will depress the ratio of cross 
sections at 630 GeV and 1.8 TeV from the pure QCD expectation. An explicit 
calculation with CTEQ4M parton densities and the b rapidity selection \y\ < 1, 
yields a pure QCD prediction at NLO of 0.17 +/- 0.02 for pffi = 10.5 GeV, and 
0.16 +/- 0.02 after inclusion of the gluino pair contribution. Either of these numbers 
is consistent with forthcoming data from CDF on this ratio.0 



2.2. Same-sign to Opposite-sign Leptons 

Since the g is a Majorana particle, its decay yields both quarks and antiquarks. 
Gluino pair production and subsequent decay to b's will generate bb and bb pairs, as 
well as the bb final states that appear in QCD production. When a gluino is highly 
relativistic, its hclicity is nearly the same as its chirality. Therefore, selection of 
g's whose transverse momentum is greater than their mass will reduce the number 
of like-sign b's. In the intermediate pt region, however, the like-sign suppression is 
reduced. The cuts chosen in current hadron collider experiments for measurement of 
the ratio of like-sign to opposite-sign muons result in primarily unpolarized g's, and, 
independent of the b mixing angle, an equal number of like-sign and opposite-sign 
b's is expected at production. The SUSY mechanism leads therefore to an increase 
of like-sign leptons in the final state after semi-leptonic decays of the b and b quarks. 
This increase could be confused with an enhanced rate of B — B° mixing. 

Time-integrated mixing analyses of lepton pairs observed at hadron colliders 
are interpreted in terms of the quantity \ = fdXd + fsXs, where fd and f s are the 
fractions of B® and B° s hadrons in the sample of semi-leptonic B decays, and Xf 
is the time- integrated mixing probability for B®. Conventional bb pair production 
determines the quantity LS C = 2%(1 — x), the fraction of bb pairs that decay into 
like-sign leptons. The SUSY mechanism leads to a new expression 

LS = \ °y + LS C J" = 2x cff (l - Xoff). (1) 

Z(Tgg + (Tq C d &gg +&qcd 

The factor 1/2 arises because N(bb + bb) ~ N(bb) in the SUSY mechanism for the 
selections on pTb made in the CDF run I analysis. Defining G = <Jgg/a qc d, the 
ratio of SUSY and QCD bottom-quark cross sections after cuts, and solving for the 
effective mixing parameter, one obtains 



x i l L J 



VTTg 2 [ VTTg 



(2) 



The CDF measurementQ of XcS = 0.131_± 0.02 ±0.016 is marginally larger than the 
world average value x — 0.118 ± 0.005jla assumed to be the contribution from the 
pure QCD component only. 
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The ratio G is determined in the region of phase space where the measurement is 
made,Q with both final 6's having prb > 6.5 GeV and rapidity \yb\ < 1. With = 
3.5 GeV, G — 0.37 and 0.28 for gluino masses nig — 14 and 16 GeV, respectively. 
The predictions are Xeff = 0.17 ± 0.02 for m g = 14 GeV, and Xcff = 0.16 ± 0.02 
with m,g — 16 GeV. Additional theoretical uncertainties arise because there is no 
fully differential NLO calculation of gluino production and subsequent decay to b's. 
The choice wig > 12 GeV leads to a calculated XcS consistent with the data within 
uncertainties. With (Jgg/a qc d ~ 1/3, the mixing data and the magnitude and px 
dependence of the b production cross section can be satisfied. 

3. Other Experimental and Theoretical Constraints 

An early study by the UAl Collaboration!!!! excludes g's in the mass range 4 < 
nig < 53 GeV, but it starts from the assumption that there is a light neutralino 
Xi whose mass is less than the mass of the gluino. The conclusion is based on the 
absence of the expected decay g —> q + q+ /Et, where /Et represents the missing 
energy associated with the Xi- I n t ne scenario discussed above, this decay process 
does not occur since the bottom squark is the LSP, the SUSY particle with lowest 
mass, and the Xi mass is presumed to be largefj.e., > 50 GeV). An analysis of 2- 
and 4-jet events by the ALEPH collaborationt-3 disfavors g's with mass nig < 6.3 
GeV but not g's in the mass range relevant for the SUSY interpretation of the 
bottom quark_production cross section. A similar analysis is reported by the OPAL 
collaboration.!^ A light b isnot excluded by the ALEPH analysis. The exclusion 
by the CLEO collaboratioiH of a b with mass 3.5 to 4.5 GeV does not apply since 
their analysis focuses only on the decays b — > civ and b — > cJL. The b need not decay 
leptonically nor into charm. On the other hand, these data might be reinterpreted 
in terms of a bound on the R-parity violating lepton-number violating decay of b 
into c£. It would be interesting to study the hadronic decays b — > cg—svith q = d 
or s, and b — > us with the CLEO data. The DELPHI collaboration'stll search for 
long-lived squarks in their 77 event sample is not sensitive to mg < 15 GeV. The 
combined ranges of b and g masses are also compatible with renormalization group 
equation constraints and the absence of color and charge breaking minima in the 
scalar potentialJlS 

There are important restrictions on couplings of the bottom squarks from precise 
measurements of Z° decays. A light b would be ruled out unless its coupling to the 
Z° is very small. The squark couplings to the Z° depend on the mixing angle 
8b- If the light bottom squark (61) is an appropriate mixture of left-handed and 
right-handed bottom squarks, its lowest-order (tree-level) coupling to the Z Q can 
be arranged to be smalLj if sin 2 8b ~ 1/6. The couplings Z^ i and Z^i survive, 
where 62 is the heavier bottom squark. However, as long as the combination of the 
masses mg + mg 2 is less than the maximum center-of-mass energy explored at LEP, 
these couplings present no difficulty. This condition implies roughly mg > 200 GeV. 
However, much lower masses of 62 might be tolerated. A careful phenomenological 
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analysis is needed of expected 62 decay signatures, along with an understanding 
of detection efficiencies and expected event rates, before one knows the admissible 
range of masses consistent with LEP data. At higher-order, unless the 62 mass 
is of order 100 GeV, contributions from loop processes in which light gluinos are 
exchanged may produce significant deviations from measurements of the ratios Ab, 
the forward-backward b ^symmetry at the Z , and Rb, the hadronic branching ratio 
of the Z° into b quarksMiJ 

Bottom squarks make a small contribution to the inclusive cross section for 
e + e~ — * hadrons, in comparison to the contributions from quark production, and 
bb* resonances are likely to be impossible to extract from backgrounds E3 The angu- 
lar distribution of hadronic jets produced in e + e~ annihilation can be examined in 
order to bound the contribution of scalar-quark production. Spin-1/2 quarks and 
spin-0 squarks emerge with different distributions, (1 ± cos 2 iD, respectively. The 
angular distribution measured by the CELLO collaborations is consistent with 
the production of a single pair of charge- 1/3 squarks along with five flavors of 
quark-antiquark pairs. Greater statistics would be valuable. 

4. Predictions and Implications 
4.1. T Decay into Bottom Squarks 

If the bottom squark mass is less than half the mass of one of the Upsilon states, 
then Upsilon decay to a pair of bottom squarks might proceed with sufficient rate 
for experimental observation or exclusion of a light bottom squark. The expected 
rate for T —* bblm&y be computed as a function of the masses of the bottom squark 
and the gluino.E 4 ] 

The data sample is largest at the T(AS). For a fixed gluino mass of 14 GeV, the 
branching fraction into a pair of bottom squarks is about 10~ 3 , for mr = 2.5 GeV, 
and about 10 -4 for m~ h = 4.85 GeV. A sample as large as 10, 000 may be available 
in current data from runs of the CLEO detector. 

The predicted decay ratescJ for the T(nS), n = 1, 3 can be interpreted as 
predictions of the width for the corresponding values of mi and rrig, or as lower 
limits on the sparticle masses given known bounds on the branching fractions. The 
current experimental uncertainties on the hadronic widths of the T's are compatible 
with the range of values of mg and rrig favored in the work on the bottom quark 
production cross section in hadron reactions described above. The analysis^ 4 ! of 
T(nS) decays shows nevertheless that tighter experimental bounds on the bottom 
squark fraction are potentially powerful for the establishment of lower bounds on 
mi and rrig. 

In conventional QCD perturbation theory, the hadronic width of the T is cal- 
culated from the three-gluon decay subprocess, T — * 3g, and T^g cx a^. The SUSY 
subprocess adds a new term to the hadronic width from T — > b + b* . If this new 
subprocess is present but ignored in the analysis of the hadronic width, the true 
value of a s (p = nib) will be smaller than that extracted from a standard QCD fit 
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by the factor (1 — rsusY/r3g)^- For a contribution from the bb* final state that is 
25% of r T (lS'), the value of a s extracted will be reduced by a factor of 0.9, at thn 
lower edge of the approximately 10% uncertainty band on the commonly quotedE3 
value of a s (mb). A thorough analysis would require the computation of NLO con- 
tributions in SUSY-QCD to both the 3g and bb* amplitudes and the appropriate 
evolution of a s (/i) with inclusion of a light gluino and a light bottom squark. 

Direct observation of Upsilon decay into bottom squarks requires an under- 
standing of the ways that bottom squarks may manifest themselves, discussed in 
more detail below. Possible baryon-number-violating R-parity-violating decays of 
the bottom squark lead to u + s; c + d; and c + s final states. These final states of 
four light quarks should be distinguishable from conventional hadronic final states 
mediated by the three-gluon intermediate state. For example, a greater rate of 
baryon antibaryon production is likely. If the b lives long enough, it will pick up a 
light quark and turn into a B-mesino, B. Charged B-mesino signatures in T decay 
include single back-to-back equal momentum tracks in the center-of-mass; measur- 
ably lower momentum than lepton pairs ( < 4 GeV/c vs. ~ 5 GeV/c for muons 
and electrons); 1 + cos 2 9 angular distribution; and ionization, time-of-flight, and 
Cherenkov signatures consistent with a particle whose mass is heavier than that 
of a proton. At stake is discovery, or new limits on the mass, of the b as well as 
measurement of or new limits on the R-parity violation couplings of the b. 

Pseudoscalar rjb decay into a pair of bottom squarks is forbidden, but the higher- 
order process of 77b decay into a pair of B-mesinos can proceed. Decays of the Xbo 
and Xb2 into a pair of bottom squarks are allowed. Implications of a lowjnass b and 
low-mass g for rare B decay phenomena are explored by Becher et alEQ 



4.2. Hadron Reactions 

Among the predictions of this SUSY scenario, the most clearcut is pair production 
of like-sign charged B mesons at hadron colliders, B + B + and B~B~ . To verify 
the underlying premise, that the cross section exceeds expectations of conventional 
perturbative QCD, a new measurement of the absolute rate for b production in 
run II of the Tevatron is important. A very precise measurement of x m run II is 
desirable. Since the fraction of Vs from gluinos changes with pxb, & change of \ is 
expected when the cut on pxb is changed. The b jet from g decay into bb will contain 
the b, implying unusual material associated with the b in some fraction of the bb 
data sample. The existence of light 6's means that they will be pair-produced in 
partonic processes, leading to a slight increase (~ 1%) in the hadronic dijet rate. 

The SUSY approach increases the b production rate at HERA and in 77 col- 
lisions at LEP by a small amount, not enough perhaps if early experimental in- 
dications in these cases are confirmedJj'Q Full NLO SUSY-QCD studies should be 
undertaken. In these two cases, the apparent discrepancy may find at least part 
of its resolution in the fact that bb production occurs very near threshold where 
fixed-order QCD calculations are not obviously reliable. Uncertain parton densities 
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of photons may play a significant role. 

4.3. Running of a s 

The presence of a light gluino and a light bottom squark slow the running of the 
strong coupling strength a s (u). Above gluino threshold, the (3 function of (SUSY) 
QCD is 

The b (color triplet scalar) contributes little to the running, equivalent to that of 
1/4'th of a new flavor, but the g (color octet fermion) is much more significant, 
equivalent to 3 new flavors of quarks. A precise determination of /3(a s ) appears 
to be the best way to confirm or exclude the possible existence of a light gluino. 
Using a method that relies heavily on a "renormalization group invariant" (RGI) 
technique to minimize non-perturbative and invers£=power contributions and scale 
dependence, members of the DELPHI collaborationEO extract the value n e ^ s = 4.75± 
0.44 from an analysis of data on the thrust distribution < 1 — T > in e + e~ — > 
hadrons over the energy range E cm = 15 to 200 GeV. A similar value for nf with 
somewhat larger uncertainties may be deduced from fitscS to the Q 2 variation of 
a s (Q 2 ). These results can be compared to the pure QCD expectation rif = 5 below 
tt threshold. The small quoted uncertainty on n e j S would preclude a light gluino. 
However, it remains crucial to understand the assumptions and constraints inherent 
in the use of the RGI approach and to verify whether the same method applied to 
other event shape- variables yields consistent results for n e ^ . 

In the standard model, a global fit to all observables provides an indirect 
measurement of a s at the scale of the Z boson mass Mz- The value a s (Mz) — 
0.1184 ± 0.006 describes most observables properly.c3 Extrapolation from Mz to 
a lower scale /i with inclusion of a light gluino reduces a s (/x) from its pure QCD 
value. The presence of a light gluino, with or without a light bottom squark, also 
requires reanalysis of the phenomenological determinations of a s (/Lt) at all scales to 
take into account SUSY processes and SUSY-QCD corrections to the amplitudes 
that describe the relevant processes. To date, a systematic study of this type has 
not been undertaken, but, as mentioned above, consistency is achieved for T decays. 
A lesser value of a s (rrib) leads, under slower evolution, to the same a s (Mz). 

4.4. b-onia 

Bound states of bottom squark pairs could be seen as J p = + , 1~, 2 + , ... mesonic 
resonances in 77 reactions and in pp formation, with masses in the 4 to 10 GeV 
range. They could show up as narrow states in the [i~ invariant mass spectra 
at hadron colliders, between the J "/\& and T. At an e + e~ collider, the intermediate 
photon requires production of a J PC = 1~~ state. Bound states of low mass squarks 
with charge 2/3 were studied with a potential model.c3 The small leptonic widths 
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were found to preclude bounds for m~ > 3 GeV. For bottom squarks with charge 
— 1/3, the situation is more difficult. 

4.5. b lifetime and observability 

Strict R-parity conservation in the MSSM forbids b decay unless there is a lighter 
supersymmetric particle. R-parity-violating and lepton-number-violating decay of 
the b into at least one lepton is disfavored by the CLEO dataEl and would imply the 
presence of an extra lepton, albeit soft, in some fraction of b jets observed at hadron 
colliders. The baryon-number- violating R-parity-violating (ft p ) term in the MSSM 
superpotential is W# p = \"j k U£DjDl; Uf and D\ are right-handed-quark singlet 
chiral superfields; and i,j,k are generation indices. The limits on individual /2 p 
and bar yor i-rmrnher violating couplings A" are relatively weak for third-generation 
squarksjifi \'/ jk < 0.5 to 1. 

The possible ft p decay channels for the b are 123 : b — > u + s; 213 : b — > c + d; 

and 223 : b — > c + s. The hadronic width isE3 

jji — 

F(6 - jet + jet) = sin 2 % £ |A|; 3 | 2 . (4) 

j<k 

If m l = 3.5 GeV, r(6 ij) = 0.08|A^ 3 | 2 GeV. Unless all A£ 3 are extremely small, 
the b will decay quickly and leave soft jets in the cone around the b. Bottom-quark 
jets containing an extra charm quark are possibly disfavored by CDF, but a detailed 
simulation is needed. 

If the b is relatively stable, the b could pick up a light u or d and become a B~ 
or B° "mesino" with J = 1/2, the superpartner of the B meson. The mass of the 
mesino would fall roughly in the range 3 to 7 GeV for the interval of b masses favored 
by the analysis of the bottom quark cross section. The charged mesino could fake a 
heavy muon if its hadronic cross section is small and if it survives passage through 
the hadron calorimeter and exits the muon chambers. Extra muon-like tracks would 
then appear in a fraction of the bb event sample, but tracks that leave some activity 
in the hadron calorimeter. The mesino has baryon number zero but acts like a 
heavy proton or antiproton - perhaps detectable with a time-of-flight apparatus. A 
long-lived b is not excluded by conventional searches at hadron and lepton colliders, 
but an analysisEil similar to that for g's should be done to verify that there are no 
additional restrictions on the allowed range of b masses and lifetimes. 
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